function [FDI] = Fdi_Calc(Amp,x_fft)
%FDI_CALC 频域特征计算
%   此处显示详细说明
% [Amp,Pha,x] = getAmpSpec(data,fs);
x_fft = x_fft';
fre_line_num = max(size(Amp));
FDI.Ampavg = mean(Amp);                                    % 频域幅值均值
FDI.Ampvar = sum((Amp-FDI.Ampavg).^2)/(fre_line_num-1);    % 幅值方差
FDI.Ampstd = sqrt(FDI.Ampvar);                             % 幅值标准差
FDI.Ampskewness = sum(Amp-FDI.Ampavg).^3/((fre_line_num-1)*FDI.Ampstd^3);    % 幅值偏度
FDI.Ampkurtosis = sum(Amp-FDI.Ampavg).^4/((fre_line_num-1)*FDI.Ampstd^4);    % 幅值峭度

FDI.fc = sum(x_fft.*Amp)/sum(Amp);             % 重心频率
FDI.Msf = sum((x_fft.^2).*Amp)/sum(Amp);       % 均方频率
FDI.rmsf = sqrt(FDI.Msf);                      % 均方根频率
FDI.vf = sum((x_fft-FDI.fc).^2.*Amp)/sum(Amp); % 频率方差
FDI.rvf = sqrt(FDI.vf);
sigma = sum((x_fft-FDI.fc).^2.*Amp)/(fre_line_num-1);
FDI.fskewness = sum((x_fft-FDI.fc).^3.*Amp)/(fre_line_num*sigma.^3);        % 频率歪度
FDI.fkurtosis = sum((x_fft-FDI.fc).^4.*Amp)/(fre_line_num*sigma.^4);        % 频率峭度
FDI.sqrtf = sum(sqrt(abs(x_fft-FDI.fc)).*Amp)/(sqrt(sigma)*fre_line_num);   % 平方根比率
end

